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We study the effects of Fermi-liquid interactions on quasi- two-dimensional d-wave superconductors 
in a magnetic field. The phase diagram of the superconducting state, including the periodic Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state in high magnetic fields, is discussed for different strengths 
of quasiparticle many-body interactions within Landau's theory of Fermi liquids. Decreasing the 
Fermi-liquid parameter Fq causes the magnetic spin susceptibility of itinerant electrons to increase, 
which in turn leads to a reduction of the FFLO phase. It is shown that a negative Fq results 
in a first-order phase transition from the normal to the uniform superconducting state in a finite 
temperature interval. Finally, we discuss the thermodynamic implications of a first-order phase 
transition for CeCoIns. 
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I. INTRODUCTION 

The coexistence of magnetism and superconductivity 
is generally thought to be mutually exclusive. Recent dis- 
coveries of exotic superconductors, which exhibit coexis- 
tence with ferromagnetism or antiferromagnetism, have 
challenged this long standing belief J^-i^iS, However, it has 
been known for some time that under certain circum- 
stances conventional superconducting order can coexist 
with paramagnetic order in high magnetic fields. Since 
the 1960s this state of coexistence has become known as 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) statc^^ More 
recently, it has been thought to be observed in the heavy- 
fermion superconductor CeCoInsi^i^ 

Here, we study the Pauli paramagnetic dcpairing effect 
of a magnetic field in a d-wave superconductor. We use 
the weak-coupling Bardeen-Cooper-Schrieffer theory of 
superconductivity with quasiparticle interactions as de- 
scribed by Landau's theory of Fermi liquids. The FFLO 
state of a spin-singlet superconductor in high magnetic 
fields appears as a result of the intricate interplay of two 
effects: (1) loss of superconducting condensation energy 
in a magnetic field when Cooper pairs (with antiparallel 
spins) are depaired, and (2) gain of magnetic energy due 
to the Zeeman effect by spin-polarizing quasiparticles. 
As a result, in high fields a spatially nonuniform super- 
conducting state coexists with pockets of spin-polarized 
quasiparticles localized at the zeros of the oscillating or- 
der parameter in real space. 

We address the modifications of the FFLO state by 
Fermi-liquid (FL) interactions. Our analysis is for quasi- 
two-dimensional (quasi-2D) fermionic systems with a 
cylindrical Fermi surface and the magnetic field applied 
perpendicular to the axis of the cylinder. In this geom- 
etry we can neglect any orbital effects on the supercon- 
ducting condensate due to magnetic field. To further sim- 
plify our analysis, we restrict our study to FFLO states 



with ID spatial modulations of the order parameter and 
neglect any low-temperature transition to a FFLO state 
with 2D modulationsi^ 

The effects of FL interactions in quasi-2D s-wave 
superconductors were first studied by Burkhardt and 
Rainer within quasiclassical theoryiiS They reported a 
considerable change of the standard FFLO phase dia- 
gram when tuning the FL parameter Fq. Since the un- 
derlying physics of Pauli dcpairing is the same for all 
spin-singlet superconductors, we expect similar new phe- 
nomena to occur for c?-wave pairing states. However, we 
anticipate additional effects due to gap nodes and the 
associated spin-polarized nodal quasiparticles. 

Here, we extend our earlier workii by including many- 
body interactions in the form of FL effects, which enter 
the quasiclassical Eilcnberger equatioufi^ii^ 

[i£mJ3 ~vz - CTfl - a, + iv/ • V g = , (1) 

with Matsubara frequency £,„ and quasiclassical Green's 
function g, through the FL dressed Zeeman term, vz, and 
the FL self-energy, (jpLr^'^ 

The Pauli matrices ai describe the coupling of quasipar- 
ticle spins to the FL dressed external field Bo/(l + Fq) 
and the internal exchange field i/. The latter satisfies the 
self-consistency condition given by the spin part of the 
diagonal component of g, 

iy{K) ^ A^oTY. f #'g(R,p';e„). (3) 

Em 

Aq is the isotropic channel of the antisymmetric part of 
the Landau interaction A"(p,p'). The Landau parame- 
ter Aq is related to the quasiparticle FL parameter Fg 
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FIG. 1: (Color online) The phase diagram of a 2D d^2_y2- 
wave superconductor for Fq = 0.5 (left panel) and Fq — 0.0 
(right panel). For positive Fq the LO state is stabilized over 
a wider range of fields. Note that the energetically unphysical 
Pauli-limited transition is first order (dot-dashed black line) 
at low temperatures, Q < T < Tp < Tfflo ~ O.SGTc. Above 
Tp ~ QATc the instability would become second order (dotted 
line). Without Fermi liquid effects Tp — Tfflo- At the 
lower critical field Bci a LO state with modulation q along 
nodal directions (110) is stablized (solid magenta line). At 
even higher fields a LO state with q || (100) becomes stable 
(dashed green line). 



through ^ F§/{\ + F^). fi = ig/2)\fiB\ is the mag- 
netic moment of an electron. The g factor of a free elec- 
tron is (7 = 2. Here g is a free parameter. 



II. RESULTS AND DISCUSSION 

In Figs. [T] and [5] we show the computed phase dia- 
grams of a 2D d-wave superconductor for tlrree different 
strengths of the FL parameter Fq ranging from negative 
to positive. The evolution of the c?-wave phase diagram 
with Fq is similar to the s-wave caseJ^ 

We determine the order of a phase transition by cal- 
culating the jump in the spin magnetization (density of 
the magnetic moment))^ 

M(R) = ^^(^^Bo-r|] y'#'g(R,p';£„)j , 

(4) 

across the transition line. Simultaneously, we check it 
by directly evaluating the free energy. A discontinuity 
of the magnetization, M = d{F/V)/d'B, defines a first- 
order phase transition, while a kink in M (discontinuity 
in the susceptibility) defines a second-order transition. 

We adopt the following notation for drawing transi- 
tion lines: second-order transitions have solid lines, while 
first-order transitions have dashed lines. For compari- 
son, the unphysical part of the normal (N) to uniform 
superconducting (USC) state transition (Pauli limited) 
is shown by a thin dot-dashed line inside the Larkin- 



FIG. 2: (Color online) The phase diagram of a 2D d-wave su- 
perconductor for Fq = -0.5 with Tfflo < Tp f» 0.75Tc. 
Note the first-order LO-N transition between Tt < T < 
Tfflo- Inset: Sketch of evolution of the order parameter 
near Bd transition from uniform (USC) to periodic (PER) 
LO solution at T/Tc — 0.15. In a narrow wedge of magnetic 
fields single (SDW) and double (DDW) domain wall solutions 
are favored over either USC or PER states. 



Ovchinnikov (LO) phase, as well as the correspond- 
ing second-order phase transition between USC and LO 
phases (solid magenta line). 

a. Second-order transition line at Bc2- First, we 
consider the second-order instability line of the upper 
critical field Bc2 from the N state into the USC state or 
from the N state into the spatially nonuniform (periodic) 
FFLO state with an order parameter A(R) ~ exp(zq-R) 
or A(R) ~ cos q • R. The phase transition can be ob- 
tained by linearizing the Eilenberger equation in A. Near 
a second order transition u is zero in linear order. Thus, 
the linearized gap equation for A is identical to that of 
the Fq = case if one replaces Bq with Bq/{\ -f- Fq ). 
So one obtains the second order normal-state instability 
line from the known solution^ i^^i^^'^^i^^ by simple scaling, 
mSc2(T; FS) = t^B,2{T; 0)(1 + Fg). 

b. First-order transition line at Bc2'- For first-order 
transitions, we must solve the general expressions ([1]) - 
([5]), which are nonlinear in the mean fields, and calcu- 
late the corresponding Green's functions and free energy. 
However, the calculation of the Pauli-limited transition 
line from the normal into the uniform state is straightfor- 
ward, since we know already the general form of g for a 
uniform superconductor j^i At T = 0, the free energy den- 
sity can be expressed in a very intuitive way, AF/V 



^AM • Bq, similar to the result by 



Clogstoufi^ where AM = Mat — M and Np is the den- 
sity of states per spin at the Fermi energy. It follows that 
a gain in condensation energy happens at the expense of 
the magnetic energy of the spin-polarized quasiparticles, 
which is proportional to the difference of the spin mag- 
netization between the N and USC state. In the s-wave 
superconductor in the absence of the Meissner effect, the 
electron magnetization M vanishes at T = 0, and the 
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Pauli limited field is fiBp ^ ^{1 + F^){\A{p- Bp)\^)^s- 
For d-wave, M(T ~ 0) is nonzero due to nodal quasiparti- 
cles, but is reduced from the normal-state magnetization 
by a fraction p = |aM|/|Mjv|. Then the right-hand side 
of the previous equation needs to be divided by ^/p. 

If many-body interactions, like FL effects, are consid- 
ered in the N state, then the spin magnetization of an 
isotropic Pauli paramagnet is given by Mjv = xn'^o ~ 
2/^^A^fBo/(1 + Fq). Thus, for positive FL parameters 
Fq the normal-state susceptibility xtv is suppressed com- 
pared to a noninteracting Fermi gas, while for negative 
Fq it is enhanced. The FFLO state, as well as the 
use state, become stable in a wider range of magnetic 
fields for positive Fq (Fig. [J), and in a smaller region 
for negative Fq (Fig. This happens because pock- 
ets of polarized electrons do not gain as much in mag- 
netic energy for Fq > as for Fq = 0, and the opposite 
happens for Fq < 0. The FFLO state exists for any 
Fq > 0, but disappears at some critical negative value, 
when the upper critical field Bc2 of the FFLO state, 
Bc2{T;Fg) = Sc2(T;0)(l drops below the Pauli- 

limited field, Bp{T;F^) w Bp{T;0)y/l + F^. At that 
point the FFLO state becomes unstable against the USC 
state for any field and temperature. The numerically 



determined critical value is Fn 



-0.765 for a d-wave 



superconductor, which is lower than that for an s-wave 
superconductor, Fq — — 0.5J^ 

We see that the FFLO state in a 2D d-wave super- 
conductor is more stable against FL effects compared to 
the s-wave case discussed in detail by Burkhardt and 
Rainerfiii Again, this is not completely unexpected, be- 
cause spin-polarized nodal quasiparticles can gain mag- 
netic energy without the breaking of Cooper pairs, which 
is unavoidable in a fully gapped superconductor. 

Further, negative values of Fq make part of the LO-N 
transition, < T < TppLO, to be first order (Fig. [2]). 
This is in qualitative agreement with the s-wave caseJ^ 
However, in contrast to the s-wave case we failed to de- 
tect any first-order transition inside the LO phase be- 
tween phases with different periods q. 

c. Second-order transition line at Bd-' Fig. [5] also 
shows details of the lower critical transition Bd from the 
USC state to the periodic LO state. We calculate the 
free energy as a function of the field for four different 
types of order parameters. We find successive transitions 
from the USC solution to the single-domain wall (SDW) 
and then to the double-domain wall (DDW) and next to 
the periodic solution (PER) in a thin but finite wedge 
of magnetic fields. For T/Tc ~ 0.15, we sketch in the 
inset of Fig. [5] the different energetically favorable solu- 
tions with their respective transitions in magnetic field. 
Although the transition from the USC to SDW state is 
continuous, there is a sequence of transitions, as domain 
walls enter the bulk one by one with increasing field. 

d. The USC- N transition: Finally, in Fig. [3] we show 
the normalized Pauli limited transition between the nor- 
mal and uniform superconducting state for several Fermi- 
liquid parameters. The break in line from solid to dashed 



indicates, as before, the change from a second to first- 
order transition. On the right side of Fig. [3] we show the 
normalized magnetization in the USC state as a function 
of field. Note that the magnetization jump, when cross- 
ing into the normal state, is larger for negative Fermi- 
liquid parameters. 
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FIG. 3: (Color online) Left panel; The Pauli-limited 
upper critical field for Fermi-liquid parameters Fq = 
{0.5,0.0,-0.5}. Right panel: Normalized uniform magnetiza- 
tion vs. field for temperatures T/T^ = {0.1,0.3,0.5,0.7,0.9} 
as labeled in the bottom window. 



III. THERMODYNAMIC IMPLICATIONS 

So far we discussed specific results for the phase dia- 
gram of a quasi-2D Fermi liquid model in the supercon- 
ducting state. As we have seen, the existence and extent 
of a first-order phase transition between the normal and 
superconducting state - USC or FFLO - at lower tem- 
peratures can be modeled by invoking a negative Fermi 
liquid parameter Fq . On the other hand, in three dimen- 
sional superconductors of type-II the first-order transi- 
tion line can be modified by a combination of Zeeman and 
orbital depairing. Generally, the existence of a first-order 
phase transition between the normal and superconduct- 
ing state is an anomalous phenomenon for strong type-II 
superconductors. Since the first-order transition is seen 
in CeCoIns, we want to address the thermodynamic con- 
straints imposed along such transition independent of a 
specific model. This is an important question that needs 
to be addressed if one wants to compare theory with ex- 
periments. 

Although most experiments report phase diagrams for 
CeCoIns in fair agreement with each other (see Refs. 
[^j8,20,21,2^, there is a noticeable variation in the po- 
sition and sharpness of the first-order transition between 
the normal and superconducting state. For that rea- 
son, we check the internal consistency of independent 
experiments by a thermodynamic analysis. Along the 
first-order transition line the free energy is continuous 
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and results in a generalized Clausius-Clapeyron equa- 
tion, = — y^j^j , which relates the jumps in entropy, 
AS* = Sn — 5*, and magnetization, AM ~ Mn — M, 
and volume V. If the magnetization in the supercon- 
ducting state is reduced by a fraction < p < 1 from 
Mat, then AS/V = -pMiydBp/dT. Consequently, the 
latent heat associated with this transition is Q = TaS = 
pTMN\dBp/dT\. 

The experiments by Bianchi et ali^ show a value of 
A^/Knoi « 200 mJ/(mol K) at T w 0.5K, and a mea- 
sured slope of the upper critical field of dBp/dT sa — 1.5 
T/K at T « 0.5 K,^i^ Whereas Tayama et al3 reports a 
magnetization jump of AM « 0.1 M^ « 80 mJ/(mol T) 
at T « 0.45 K. It is obvious that the agreement between 
the ratio of the discontinuities, AS/{VaM) ~ 2.5 T/K, 
and the slope of Bc2 is poor at T w 0.5 K (a deviation 
of ~ 70%), despite good overall agreement between both 
phase transition lines. The origin of this inconsistency is 
poorly understood, but might be related to the nature of 
the localized / electrons and their contributions to the 
magnetization. Further entropy and magnetization stud- 
ies are needed to resolve this open problem. 



IV. CONCLUSIONS 

We studied the superconducting phase diagram of 
quasi-2D d-wave superconductors in the presence of 
Fermi- liquid effects in high magnetic fields. We found 
that a negative Fermi-liquid parameter Fq increases the 
gain in magnetic energy, while at the same time it reduces 
the available phase space of a FFLO state. The uniform 
superconducting and periodic FFLO state are competing 
and at a critical Fermi-liquid parameter Fq « —0.765 the 
FFLO state is completely suppressed. 



We note that in order to explain the high-field phase 
diagram of the CeCoIns superconductor in terms of an 
FFLO state one needs to go beyond the simplistic Fermi- 
liquid picture considered here. While we find that the 
inclusion of Fermi-liquid effects considerably changes the 
phase diagram, the changes are not consistent with the 
experimental findings. For example, (1) the magnitude of 
the calculated critical temperature Tp, where the Pauli- 
limited upper critical transition changes from second to 
first order between the uniform superconducting and nor- 
mal state, and (2) the corresponding magnetization jump 
are much larger than those seen in experiment. We also 
note that the shape of the transition lines of the calcu- 
lated FFLO state is qualitatively different from exper- 
iments. It is not unreasonable to expect that some of 
those discrepancies may be overcome by including the 
effects of impurity scattering,— antiferromagnetic spin 
fluctuations, local magnetic moments, or orbital effects. 

Finally, our analysis of the first-order transition puts 
stringent constraints on thermodynamic properties and 
reveals a significant discrepancy between specific heat 
and magnetization measurements that requires further 
studies. 
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